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ABSTRACT 

It is generally considered, as a rule of thumb, that carbon monoxide forms very 
early in envelopes of AGB stars, and that it consumes most of the carbon, 
or most of the oxygen, depending on whether the photosphere is oxygen-rich 
or carbon-rich, respectively. This work focuses on the latter case , with the 
purpose of quantifying the remaining fraction of gaseous carbon which is then 
available for forming carbonaceous grains. Since AGB stars are (probably) the 
main providers of cosmic carbon grains, this residual fraction is essential in 
establishing the validity of current grain models. Here, we use a kinetic treat- 
ment to follow the chemical evolution of circumstellar shells towards steady 
state. It is shown that the residual fraction depends essentially on the atomic 
ratio of pristine gaseous carbon and oxygen, and on the cross-section for CH 
formation by collision of C and H atoms. It lies between 55 and 144 C atoms 
per 10^ H atoms, depending on the values adopted for unknown reaction rates 
and cosmic C abundance. This is much larger than predicted by the rule of 
thumb recalled above. 

The present results depend strongly on the rate of the reaction C-I-II— +CH: 
far from thermodynamic equilibrium (which is the case here), CO cannot be 
formed if this rate is as low as generally assumed. We have, therefore, es- 
timated this rate by chemically modelling the reaction and found it indeed 
much higher, and high enough to yield CO abundances compatible with ob- 
servations. An accurate experimental rate determination is highly desirable. 

Key words: Stars:Astrochemistry; circumstellar matter; AGB and post- 
AGB; stars: carbon; ISM: dust; molecular processes. 
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1 INTRODUCTION 

The spectral analysis of the CS (CircumStellar) envelopes of AGB stars has been an ac- 
tive field of observational astronomy for decades (see, for instance, Bujarrabal et al. (119941) . 
Woods et al. (120030 ). As a result, the list of detected atoms, ions, radicals and molecules is 
continuously enriched and their characterization improves in extent and quality as instru- 
ments and telescopes are constantly upgraded (see van Dishoek (120040 for references). 

One of the most abundant molecules in both 0-rich and C-rich AGBs is carbon monoxide, 
CO, so much so that the dominant isotopic variety, ^^CO, is optically thick in its character- 
istic infrared bands, precluding direct determination of its column density towards the star. 
One way of circumventing this difficulty is to measure instead the density of the rarer iso- 
tope, ^'^CO, from which that of ^^CO is deduced, assuming their ratio is somewhere between 
1/50 and 1/100, because of (undetermined) fractionation processes which locally alter the 
relative cosmic abundance of 1/89 (see Burgh at al. (120071) . Liszt (I2007p ). 

In spite of these achievements and because of observational uncertainties (such as the 
CO isotopic ratio) and the relative scarcity of stars amenable to in-depth analysis, it is 
still difficult to draw a complete and accurate observational picture of the variations of CO 
abundance with the type of AGB star. 

For want of a better guess, one may assume an arbitrary but reasonable relative abun- 
dance, such as [C0]/[if2] = 10~^ (e.g. Woods (I2003P )). Another simple rule of thumb has 
often been used: all the oxygen (respectively carbon) expelled by C-rich (respectively O- 
rich) stars is locked into gaseous CO. This rule is strongly suggested by the high reciprocal 
chemical affinities of the two elements, and the strength of the CO bond. 

Precisely because of the latter, the carbon and oxygen atoms locked in CO are irreversibly 
lost for carbonacious and silicate grains which ultimately condense in the envelope. If strictly 
applied, the rule above implies that only stars with C/0>1 can contribute to carbon grains. 
Now, for many carbon grain models, the C budget is already quite tight, so the available 
C input deserves careful scrutiny. In particular, the carrier of the extinction bump at 2175 
A (see, for instance. Bless and Savage (119720 ) requires large quantities of pure carbon (see 
Snow and Witt (I1995P for a discussion). Since CS shells, together with supernovae and novae, 
where grain condensation probably occurs in similar ways, are the main, if not sole, genitors 
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of grains, accurate and generally applicable procedures are required to estimate the amount 
of C available for grain condensation. 

On the theoretical side, a parallel effort was made to model the chemistry which takes 
place in these envelopes, using long chains of known processes and reactions. Assuming 
local thermodynamical equilibrium (LTE) has been a powerful approach for estimating the 
relative abundances of the various species of interest as a function of temperature, pressure 
and distance from the central star (see Tsuji (11964p . Salpeter (11974^ . Duley and Williams 
(I1984P ). However, for LTE to be applicable, steady states at high densities and temperatures 
are required, so that detailed balance may be invoked, which is not the case in CS. Besides, 
it is not easy to include solid grains in the LTE scheme. 

Another approach consists in selecting a number of molecular species and a number of 
known reactions between these species, and expressing the rate of change of each species as 
a function of the instantaneous densities of all species and the rates of their reactions. In 
general, a steady state is quickly established, which depends on local temperature and initial 
densities. The complexity of the problem is illustrated by the work of Prasad and Huntress 
( I1980p . a milestone in the study of the chemistry of intersellar molecular clouds. While the 
physical conditions in clouds are not comparable with those of CS envelopes, that work is 
an inspiring example for our present purposes. Cherchneff ( 12006P recently used the kinetic 
approach to study the dependance of molecular populations in the inner parts of CS of AGB 
stars, as a function of the C/0 ratio in the photosphere. The present research is also based 
on chemical kinetics rather than LTE, but the goal is much less ambitious, being limited to 
estimating the fraction of carbon available for grain formation. 

We consider carbon grains and gaseous CO to be competitors for the available C atoms 
from the photosphere. A small number of chemical reactions between the most abundant 
gaseous chemical species are retained to define the two corresponding routes for pristine 
C atoms (Sec. 2). They translate into a set of differential equations which can be solved 
numerically for given relative abundances of C and O, to yield the relative number of C 
atoms that end up in grains and CO respectively (Sec. 3). The results are arranged as 
a function of the ratio C/0, displaying continuous trends from 0-rich to C-rich stars. The 
fraction of carbon available for grains is thus predicted to be overall much larger than present 
estimates (Sec. 4). 

The discussion of these results (Sec. 5) singles out the reactions which dominate the 
production of CO or carbonaceous grains, and should therefore be more fully documented 
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by laboratory measurements. Finally, we discuss how the average ratio of carbon atoms 
available for grains to total H atoms can be predicted based on the present computations. 

2 THE CHEMICAL REACTIONS 

Although its bonding energy is particularly strong (250 kcal/mol or about 11 eV), CO is not 
known to form directly by simple association of neutral C and O atoms. This is analogous to 
the case of H2: because of the similarity of the two interacting atoms, conservation of energy 
and momentum requires the presence of a third body (grains or high pressure gas). This is 
not available in the environments considered here. The emission of a molecular resonance 
photon of 11 eV could help, but no such resonance is tabulated (see Radzig and Smirnov 
(I1985P ). We therefore look for indirect paths, involving successive binary reactions. 

The number of possible reactions in the CS gas is huge (cf Prasad and Huntress (I1980p . 
UMIST Data Base fl2000p ). even for the small number of atomic species that are involved 
initially: H, C and O. For present purposes, it is necessary and sufficient to select a few 
of them that could be considered as essential or sufficiently representative of the complete 
process. We now set out to do this. 

The gas densities and temperatures to be considered are characteristic of inner regions of 
AGB envelopes illuminated by photospheres at 4000 K or less, and shielded from the ionizing 
UV photons of the IS medium. Neutral radicals and molecules will therefore be dominant 
over ions, as suggested by thermodynamic equilibrium models (Tsuji ( I1964p . Duley and 
Williams (I1984p ). Of course, cosmic rays pervade CS envelopes, but their rate of ionization 
of matter (e.g. 10^^'' s^^ per H atom) is negligible, considering the rates of chemical reactions 
to be retained below (in the order of 10~^ s^^). 

Of the most abundant elements, helium will be disregarded because it is chemically 
inactive, as will be nitrogen because it does not enter in grains in quantities large enough 
to compete with the processes considered here. 

The initial reservoir of the atomic species at the start of calculations will therefore consist 
of H, C and O. Their absolute densities only affect the time constants of the reactions, not 
their outcome, which is determined by their relative densities. The H density will be set at 
10^*^ cm^^ for all computed cases, which corresponds to a mass loss rate of about 3 10~^ 
M0/y, for a wind velocity of 10 km.s~^, at a radial distance of 10^^ cm from the center 
of the star. Given the very low relative abundance of the other elements, this H density is 
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not notably altered by the chemical evolution of the gas, and may therefore be assumed to 
remain constant, for the sake of simplicity. 

In selecting possible reactants, preference will be given to those carrying smaller number 
of atoms of higher cosmic abundance. O2 is excluded on the grounds that, in this case, the 
balance between formation and destruction processes is highly unfavourable to its presence 
in significant quantities. 

Kinetic and LTE calculations as well as observations show that the abundance of gaseous 
molecules with more than 3 or 4 atoms is much lower than that of CO (see Duley and 
Williams (119841) ). Besides, the products of chemical reactions between carbonaceous reac- 
tants with more than 3 atoms are either light carbon-poor molecules (mainly CO), or heavier 
hydrocarbon molecules. This is because C can form rings or long chains, while O cannot. We 
therefore limit our chemical scheme to reactants with less than 4 atoms each, and assume 
that products with 2 or more C atoms will further react to end up mainly in carbon grains. 
This assumption is borne out by our calculations below. For instance, in Fig. 2, the only 
carbon-bearing molecules that survive to any significant extent are CO and C2H2. Reaction 
between O and C2H2 may still yield CO, but its rate is negligible (Rate06, No 321). Other 
pathways from C2H2 have been considered by Cherchneff et al. ( (I1992p : e.g. to C4H4 through 
the vinylidene radical C2HH), however, they found that, ultimately, C2H2 remained the most 
abundant molecule from which carbon grains are grown. 

Reactions with activation energy higher than 2500K are excluded, as are those with a 
rate prefactor a< 10^^^ cm^s~^. In further limiting the number of reactions to be retained 
among the selected species, we require that one at least of the products also belongs to this 
selection. Finally, trimolecular reactions are obviously excluded at the relevant gas densities. 

These (arbitrary but argumented) conditions limit the number of selected species to 11: 
H, H2, C, C2, O, OH, CO, CH, CH2, C2H and C2H2. The reactions involved are hsted in 
Table 1 with their particulars as given in the UMIST Data Base, RateOG (120071) . Interestingly, 
no way was found to form CO or carbon grains without CH as an intermediary, and the 
only route to CH, according to the cited data bases, appears to be radiative association 
of C and H. But, even with prior ionization by cosmic rays, the rate given in Rate06 is 
extremely low (10^^^ cm'^s"^). True enough, the rates of reactions 11 (H+C2 ^CH+C) and 
47 (H2+C^CH+H) are quite high, but they were excluded because of excessive activation 
energy (30450 and 11700 K, respectively). We were therefore led to use chemical modelling 
in order to study the reaction C+H— i-CH (details in the Appendix). It was found that it 
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Table 1. Reaction scheme: k(cm-'s-i) = 10-"' a(T(K)/300)''cxp(-c/T(K)); for kl, k2, see Appendix. 



RateOe No 


React ants 


Products 


a 


b 


c (K) 
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H+CH 


C+H2 


1.3 





80 


2 


H+CH2 


CH+H2 


0.664 
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H2+CH 


CH2+H 


5.46 





1943 
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C+CH 


C2+H 


0.659 
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C+CH2 


C2H+H 


1 
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C+OH 


CO+H 


1 
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C+C2H 


C3+H 


1 








138 


CH+O 


OH+C 


0.254 





2381 


140 


CH+O 


CO+H 


1.02 





914 


237 


CH2+CH2 


C2H2+2H 


1.8 





400 


241 


CH2+0 


CO+2H 


1.33 








242 


CH2+0 


CO+H2 


0.8 








316 


0+C02 


CO+C 


6 








317 


0+C2H 


CO+CH 


0.17 








kl 


H+C 


CH 


1 








k2 


CH+CH 


C2H2 


1 









has no activation energy, and the cross-section is of order 30 A ^ in the temperature range of 
interest. The same procedure was apphed to the association of CH with CH, although this 
is not as crucial as the former reaction. 

3 THE EVOLUTIONARY EQUATIONS 

The calculations will show that the transient between initial and steady states is so short 
that it takes place entirely in a very small radial span. This spatial evolution can therefore 
be overlooked. Let nniO), nc{0) and no(0) be the initial densities, and n the time derivative 
of n. Based on the essential constancy of the gas temperature and the small dependance of 
the rate constants on temperature, we shall replace each rate constant, k, by its pref actor 
written as a (cm^s~^), followed by its ranking number in RateOG. 

The following equations were solved as difference equations, with a time step dt=10~^s. 



hH2 = -a48nH2ncH + aiUHUcH + a2nHncH2 + a242noncH2; 
nc = -aesUcncH - aj^ncnoH - h iicUh 
-aG5ncncH2 + ainnncH; 
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no = -ai^oUoUcH - aissnoncH - a3nnonc2H - a^i&nonc2 
—a24,\noncH2 — o,2i2'noncH2] 

nco = a-r^ncnoH + auoUcHno + a2iincH2no + (i242'ncH2'no 
+a3i6nonc2; 

hcH = -ainnncH - amncncH - aissncnno - a-iAoncHno 
-2k2ncHncH - a4snH2ncH + a2nHncH2 + a3i7nc2Hno + kiUcUH] 
rioH = -a^^noHric + ai^sriCHno] 
nc2 = -aziQnc2no + aQsUcncH] 

ncH2 = -a2ncH2nH - amncncH2 - '2a237ncH2ncH2 - a2iinoncH2 

a2A2noncH2 + a48nH2ncH] 

nc2H = -a3i7nc2Hno - auncnc2H + amncncH2\ 
'nC2H2 — 0,237^0 H2'nCH2 + k2ncH'ncH', 
riH = nHiO) - ncH - nc2H - 2ncH2 - 2nc2H2; 
nc = nc{0) - ucH - 2nc2ff - ncH2 - 2nc2i?2; 
no = no{^) - nco - non] 

Ucgr = 2(nc2 + nc2H + nc2H2)- 



Note that, in order to ensure conservation of atoms, current values of first generation 
atoms, uq and uq-, arc not computed by incrementation like the reaction products. The 
spatial density, ncgr-i of C atoms trapped in carbonaceous grains is defined as a sum over 
C-enriched, later generation, species. The equation for uh is only meant to check that it 
does not vary significantly. 



4 THE RESULTS 

Several runs were made to cover a significant set of parameters. As an illustration. Fig. 
1 shows the time evolution of the densities of H, H2, C and O for initial values nc(0) = 
?^o(0) = 3 10^ and 77,^(0) = 10^° cm~^. Steady state is essentially reached in about one day, 
justifying the neglect of the space variable. Molecular hydrogen also forms quickly but does 
not exceed ~ 0.1 n//. 

Figures 2 and 3 show that carbon is entirely consumed into CO and grains while oxygen 
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10 10 
t (sec) 



Figure 1. A typical example of computational results. The parameters are the initial values of the spatial densities of first 



generation species (cm '^): nu 



10^ 



no 



3 10 , n//2 = 0. Although initially absent, molecular hydrogen evolves 



quickly, then levels off at a small density relative to H. Carbon is totally consumed, contrary to oxygen. Steady state is 
essentially reached in a day's time. The hydrogen density decreases only slightly, due to the formation of hydrogenated species. 



ends up partly in OH and mainly in CO. Figure 2 illustrates the detailed budget of carbon: 
C2 and the primary hydrocarbons, CH, CH2, C2H quickly decay into C2H2, which, in this 
truncated scheme of reactions, is a final product, as is CO. Now, C2H2 is known to be the seed 
of carbonaceous grains (see Frenklach and Fiegelson (I1989p ). This justifies the truncation of 
the reaction chain for our present purposes. 

Several runs with different initial densities of mother atoms showed that asymptotic 
values of products, expressed as fractions of ^c(O), vary notably with the ratio nc{0)/no{0), 
but are quite insensitive to the absolute values of these two densities. The asymptotic values 
of interest, nco{f)/nc{0) and ncgr{f)/nc{0), are therefore plotted in Fig. 4 and 5 as a 
function of nc{0)/no{0), for a few cases spanning the relevant range (filled squares; straight 
lines are drawn between squares to help the eye). 

It may help to know that a square polynomial fit to the points in Fig. 4 gives, to a very 
good approximation. 
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10S 




t (sec) 

Figure 2. Same initial densities as in Fig. 1. Evolution of C-rich species. They all decay to zero, except acetylene. CH is the 
indispensable starter and mediator of the evolution. 



> = 1.018 - 0.328 - 0.033 . 1 

In Fig. 5, the filled square with abscissa 4 is lower than the corresponding thumb- rule 
value. This is abnormal and due to insufficient time length of the run: the curve ncgr=^{^) 
did not reach its asymptote in this case. 

Finally, if one is interested in the fraction of oxygen available for silicate grains, an upper 
limit for it is no{f)/no{0), which is plotted in Fig. 6. 

5 DISCUSSION 

Back in Fig. 4 and 5, dotted lines have been drawn to illustrate the application of the rule of 
thumb recalled in the Introduction. As expected, the latter approximation is seen to agree 
with our more elaborate calculations only in the two opposite limits of nc{0)/no{0). The 
difference is considerable in between, which is the most often encountered range in the sky 
(see for instance Kwok et al. (I1997P ). 
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t (sec) 

Figure 3. Same initial densities as in Fig. 1. Evolution of O-ricli species. Oxygen is essentially consumed into CO, and partly 
OH. 

In a further effort to gauge the relevance of our chemical scheme, we sought to reduce 
the number of intermediate species to the strict minimum. Only CH was found to be in- 
dispensable for the formation of both CO and grains. The corresponding reduced scheme 
is 

he = -aesUcncH - ki ucUh + ainnncH] 

ho = -auoUoncH] 

hcgr = '2k2ncHncH + kesUcHnc; 

hcH = -ainnncH - a-mncncH - aumcHno 

-2k2ncHncH + hncnH] 

nco = no{Q) - no- 

The results are represented in Fig. 4 and 5 by open circles. These results are clearly 
incorrect considering their behaviour at the two ends of the X-axis. However, one notes that 
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Figure 4. The fraction of initial carbon tiiat is consumed in CO, as a function of the initial ratio of C to O (C/O=nc(0)/nQ(0); 
f designates the end of the calculation run). Filled squares: solutions of the set of equations in Sec. 3 (straight lines were drawn 
between computed points as an aid to the eye). Dashed line: illustration of the usual rule of thumb (see text): it runs through 
the computed points in the limits of very small and very large C/O, but the deviations are considerable in the intermediate 
range. Open circles: solutions of the reduced set of reactions and equations (Sec. 5) ; inspite of the drastic simplification, the 
general trend is similar to that of the squares, with differences smaller than a factor 2. 

^co(/)/^c(0) computed with this scheme is never less than one half of the value delivered 
by the more elaborate scheme, and the strings of squares and circles are closely similar. This 
is witness to the robustness of the first scheme. 

This was further tested by performing other runs with the first scheme, only changing one 
or two parameters. Figure 7 plots nco{f)/f^c{^) (open squares) and ncgr{f)/nc{0) (filled 
squares), for some of these cases, all with nc(0) = no(0) = 3 10^ cm~^. The integer abscissa 
refers to the sentence in the legend which describes the corresponding changes applied to the 
parameters. The results are arranged so as to higlight the opposite effects of various changes 
in the chemical scheme, with respect to the "standard" case, corresponding to the equations 
of Sec. 3. Some elements of the scheme have no great impact on the final values of interest 
here; for instance, the inclusion of OH (N=3) and the initial fraction of molecular hydrogen 
(N=5). By contrast, the association of C and H to form CH (reaction kl) is essential for 
a high yield of carbon grains, as stated above and as evidenced by the steep reduction 
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Figure 5. The fraction of initial carbon that is ultimately consumed into carbon grains (represented here by C2H2 essentially). 
Same symbol conventions as in Fig. 4. 

in carbon grain formation when kl is reduced to 0.1 10~^° (N=l). Other parameters have 
moderate impact. 

Note that, in case N=l, the yield of CO has not diminished, and even increased some- 
what, although CH is an important factor of CO production through reaction 140 (Table 
1). This is because other quite efficient, reactions are still at play, which compensate the 
blocking of that route. However, the steady state is retarded by a factor 10. 

Other computations also showed that changing all initial atom densities in the same 
proportions by a factor 2 had only marginal effects on the final yields, but retarded or 
advanced the setting of steady state. 

Finally, let us consider the problem of carbon grain formation. The fraction of cosmic 
carbon that is available for grain formation is usually equated to the difference between 
photospheric (cosmic) and ISM abundances of carbon, divided by the former. In our treat- 
ment, this is represented by the average of ncgr{f)/'nc{0) over all stars. This quantity may 
be deduced, with the help of Fig. 5, from the measurement of the photospheric C/0 ratio of 



On the formation and abundance of CO in envelopes of AGB stars 13 




Figure 6. An upper limit to the fraction of oxygen available for silicate grains, no{f)/n,o{0). 

a large number of stars of different types. // the stars were evenly distributed over the range 
0.5<C/O<2, the fraction would be ~0.25. If, furthermore, the cosmic C/H ratio is assumed 
to be that recommended by Snow and Witt fll995p . 2.25 10~^, then, the available number 
of C atoms for carbon grains becomes 55 for 10^ H atoms. 

However, this must be considered a lower limit, as it is the result of calculations with 
kl=l 10^^° cm^s~^, a reasonable but arbitrary value, adopted in Table 1, in the absence 
of experimental data. If, on the other hand, one takes kl=6.6 10^^°, as suggested by our 
chemical modeling (see Appendix), and if, moreover, one adopts the Grevesse cosmic carbon 
abundance for the Sun fll99ip . i.e. 4 10~^, then the number of available C atoms rises to 
144. 

Circumstances may be envisioned where this upper limit may still be raised to some 
extent. For instance, as the expanding gas reaches the outskirts of the envelope, it is no longer 
shielded from the UV radiation of the ISM. If the latter is strong enough, CO may then be 
dissociated, making some more free C atoms available for grain formation, as observed by 
Herpin et al. (I2002p in the case of post-AGB stars. 
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Figure 7. Effects of changes in the parameters, one at a time, defined by the rank, N. N=f: l<f=0.f 10 N=2: k2=0; 
N=3: length of run multiplied by 2. N=4: noH = at all times. N=5: all rates as in Table 1. N=6: nniiO) /"n-H = 0.1. N=7: 
kl=2 10^10. N=8: kl=6.6 10~10.N=9: a^ia = 0. N=10: 0241 +0242 = 0. C/0=1 for all. 

Another instance is a high initial ratio of molecular to atomic hydrogen. Following a 
suggestion by the reviewer of this paper, our program was run (after minor changes) for 
the case where nniO) = nH2{0) = 10^°cm~^ and nc{0) = no{0) = 3 10^ cm~^. Comparing 
with Fig. 1 to 3, the final value of ricgr was found to increase to 1.2 10^ from 0.75 10^ cm~^, 
and the final value of no to decrease to 1.8 10^ from 2.24 10^ cm~^. This notable increase 
in grain abundance in the presence of a high density of H2 is due to reactions 48 and 237 
which generate CH2 and C2H2, the progenitors of grains. 

The range 55 to 144 for 10^ H atoms, deduced from our model, is consistent with spectro- 
scopic observations of a couple of stars as reported by Snow and Witt (I1995p . in a compilation 
of measurements. While more astronomical observations could help, the present work should 
be an encouragement to further study simple chemical reactions between neutrals, especially 
reaction kl, so as to restrict the range of uncertainty of the model. Such reactions have be- 
come more amenable to measurement with the development of specialized apparata, e.g. the 
CRESU supersonic-fiow machine (see for instance Rowe (120001) ). 
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6 APPENDIX: CHEMICAL MODELING OF BINARY REACTIONS 

The rate listed in the UMIST Data Base, RateOG, for the radiative association of C and H 
to give CH is 10"^"^ cm'^s"^. Such a low value can hardly start a chain of reactions leading 
ultimately to the formation of enough carbonaceous grains to fit observations. For want of 
other experimental data, we resorted to chemical simulation. This consists in using molecular 
dynamics to represent the collision of C and H atoms, for different relative velocities and 
impact parameters, to yield the cross-section for association as a function of relative velocity. 
At the end of a run, either the two atoms fly away in opposite directions, or the CH molecule 
is formed, with internal and external energies such that the initial energy and momentum are 
conserved. An order of magnitude cross-section is then deduced from these results; combined 
with a reference temperature (and hence a velocity of impact), this yields the rate, k. 

A description of the software package used for simulation can be found in Papoular 
( I200ip . The particular code used here is AMI (Austin Model 1), as proposed, and optimized 
for carbon-rich molecules, by Dewar et al. (119851) . This semi-empirical method combines a 
rigorous quantum-mechanical formalism with empirical parameters obtained from compari- 
son with experimental results. It computes approximate solutions of Schroedinger's equation 
for each instantaneous spatial distribution of atomic nuclei, and uses experimental data only 
when the Q.M. calculations are too difficult or too lengthy. This makes it more accurate than 
poor ab initio methods, and faster than any of the latter. 

The molecular dynamics relies on the Born-Oppenheimer approximation to determine 
the motions of atoms under nuclear and electronic forces due to their environment. At every 
step, all the system parameters are memorized as snapshots so that, after completion of the 
run, a movie of the reactions can be viewed on the screen, and the trajectory of any atom 
followed all along. This makes for a better understanding of the details of mechanisms and 
outcomes. In particular, it becomes possible to measure the activation energy, if any. 

Note that the dynamics of bond dissociations and formation can only be simulated by 
using Unrestricted Hartree-Fock (UHF) wave functions in the Q.M. part of the calculation 
(see Szabo and Ostlund (I1989P ). The elementary calculation step was set at 10^^^ s. 

Semi-empirical simulation methods do not account for purely quantum-mechanical phe- 
nomena like tunneling or zero-point energy. In the present context of relatively high temper- 
atures, the former is irrelevant, and the latter hardly affects the probability and energetics of 
the physico-chemical reactions involved. The simulation does not treat absorption or emis- 
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sion of light; so the contribution of hne radiation to the energy and momentum budgets is 
not included, which is a weakness of the procedure. However, CH, for instance, has only 3 
resonance lines, all of order 3 to 4 eV (Radzig and Smirnov fll985p ). not enough to play a 
significant role in the reaction budget. 

The computations were performed for impact parameters p=0, 1, 2, 3, 4, 5 A and H 
velocities, V, between 5 and 110 A per picosecond (5 10^ and 1.1 10^ cm.s^^). No association 
was observed for p>3 A , nor for V higher than 5 10^ cm.s^^ (which corresponds to a H gas 
temperature of ~1200 K). For p=3 A , association occured for all velocities lower than 
5 lO^cm.s"^. No activation energy was detectable. 

If the cross-section radius is taken as 3 A , the area becomes 27.2 10~^^ cm^. At a gas 
temperature of 300 K (reference temperature for rates in RateOG), the relative velocity of C 
and H is 2.4 10^ cm.s~^ and kl=6.6 10^^*^ cm^s^^. However, caution is in order here, because 
the dependence on temperature of neutral-neutral reaction rates in general is not yet fully 
documented experimentally (see Smith et al. (120041) ): so kl was set at lO"^'' cm'^s"^ in Table 
1, but computations were also made for 0.1 and 6.6 10~^° cm'^s"^, for comparison. Similar 
computations were made for CH-|-CH^C2H2. 
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